Ergodic versus nonergodic behavior in oxygen deficient high- Tc superconductors 
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The oxygen defects induced phase transition from nonergodic to ergodic state in superconduc- 
tors with intragrain granularity is considered within the superconductive glass model. The model 
predictions are found to be in a qualitative agreement with some experimental observations in de- 
oxygenated high-Tc single crystals. 
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I. INTRODUCTION 

According to the recent findings (see, e.g., [1-8] and 
references therein), high-Tc superconductors (HTS) ex- 
hibit an anomalous (nonclassical [8]) magnetic-field be- 
havior, which has been attributed to the "field-induced 
intragrain granularity" in oxygen deficient samples and 
interpreted in terms of the field-induced decoupling of 
regions of oxygen-rich material by boundaries of oxygen- 
poor material. A "phase diagram" Hm{S, T), that demar- 
cates the multigrain onset as a function of temperature 
and oxygen deficiency,(S, was found [3] to confirm that 
oxygen-deficient single crystals exhibit behavior charac- 
teristic of homogeneous superconductors for H < H„i 
and inhomogeneous superconductors for H > Hm- The 
granular behavior for H > H„i has been related to the 
clusters of oxygen defects (within the CuO plane) that re- 
strict supercurrent flow and allow excess flux to enter the 
crystal. The observed Hm {S, T) data were described by a 
2-D percolation model for oxygen defects. It means that 
there exists a critical oxygen deficiency, (5c, above which 
there are no continuous current paths. For 6 greater than 
(5c, oxygen- rich superconducting "grains" are separated 
by oxygen-poor insulating boundaries so that there is 
no superconducting path through a sample. For 5 less 
than (5c, a complete current path spans the sample and 
resistance measurements show metallic behavior with a 
superconducting transition. Since Hm{S,T) signals the 
onset of granularity, a sample with i?m(<5, T) = implies 
that the crystal has so many oxygen defects that it never 
exhibits single-grain behavior [3] . 

The aim of the present paper is to show how the lack 
of oxygen in HTS materials inspires the phase transition 
from nonergodic (in nearly full-oxygenated crystals) to 
ergodic (in highly oxygen-depleted crystals) state within 
the so-called superconductive glass (SG) model (see, e.g., 
[9-15] and references therein), and to face the model pre- 
dictions with some experimental data for deoxygenated 
HTS single crystals. More specifically, the nonergodic 
(phase- coherent) state is attributed to nontrivial equilib- 
rium (long-time) behavior of the defect-free crystal (with 
(5 ~ 0) while the ergodic (paracoherent) phase corre- 
sponds to the equilibrium state of the defected crystal 
(when 5 ~ 6c). 



II. THE MODEL 

The SG model is based on the well-known Hamilto- 
nian of a granular superconductor which in the so-called 
pseudospin representation has the form [9-15] 
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where 



no = -J2JiS,T)coscj},,{H) 



(,5, T, H) = J{S, T) exp[ay [H]] , 
(^y(i?) = (^,-(^,-Ay-(i?), 
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This model describes the infinite-range interaction be- 
tween oxygen-rich superconducting grains [with phase 
4>i{t) or Josephson pseudospins Sf = e+"'^'], arranged 
in a random two-dimensional (2D) lattice (modeling the 
CuO plane of oxygen-depleted YBa2Cu307-s, where a 
glass-like picture is established [1-7]) with coordinates 
T^i = {xi,yi,0). The grains are separated by oxygen- 
poor insulating boundaries producing Josephson cou- 
pling with energy J{S,T). The system is under the in- 
fluence of a frustrating applied magnetic field iJ, which 
is assumed to be normal to the CuO plane of HTS. The 
increase of the oxygen deficiency, (5, leads to the decrease 
of the Josephson energy (via the increase of the insu- 
lating layer between oxygen-rich grains). For small S 
(such that (5 <^ 1) we can approximate the S depen- 
dence of the Josephson energy by a linear law [4] , namely 
J((5,T) w J(0,T)(1 - (5). The superconducting current 
through the Josephson junction between grains i and j 



^ 2p 7 
4(i?) = — sin(/),,(i?) 

= ^Im{J^jS+S]-} 



(4) 
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induces a diamagnetic moment of the weak-link network 
[9-13] 



N 



(5) 



To study dynamic (relaxation) behavior of the model 
(1), let us assume that in addition to the constant frus- 
trating field H, the superconducting grains are under 
the influence of a small time- varying field H\{t) <C H, 
so that cos((/)i — — Aij{H + Hi{t))) = cos (pijiH) + 
Aij{Hi{t)) sin (j)ij{H). In view of Eqs.(l)-(5), the total 
(perturbed) Hamiltonian can be cast into the form 



n{t)=Ho{H)-0iit) 



(6) 



If the perturbation is applied continuously from t = —oo 
up to f = and is cut off at t = 0, then the linear (with 
respect to the small perturbation field Hi{t) = Hi9{—t)) 
response M{t) = < fiz >/V will relax to its equilibrium 
value Meg = limj^oo M{t) according to the formula [17] 



r 

M{t) -Meq= J dt'G{t - t')Hi{t') 



(7) 



= / dt'G{t')Hi 



According to the fluctuation-dissipation theorem [17], the 
response function G(t) is related to the relaxation func- 
tion $(f) as follows, G{t) = -{d/dt)<P{t). Thus the above 
equation reads 



M{t)-M,, = -H,[m-Hoo)], (8) 



where 



$(t) =/3< M.(i)M^(0) >. 



(9) 



Here (3 — l/fcsT, the bar denotes the configurational av- 
eraging over the randomly distributed grain coordinates 
(see Appendix A) , < ... > means the thermodynamic 
averaging with the Hamiltonian 'Ho{H), and we have as- 
sumed that H = (0,0, i?) and Hi = (0,0, i?i). There- 
fore, the function <I>(t) describes the relaxation of mag- 
netization M{t) after removal of the outer disturbance. 
As a result of configurational averaging, the relaxation of 
magnetization can be approximated by the formula (see 
Appendix A) 



Mit) 
M{S,T,H,Hi) 



MiS,T,H,Hi)\D{t)\' 
x{5,T,H)Hi, 
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Here Hq = 0o/s is a characteristic Josephson field with 
s = TTC?^ an average JJ projection area, and N is the num- 
ber of grains. Thus, all information about the dynamic 
(relaxation) properties of the system is contained in 
the time-dependent correlator D{i) = Dij{t) 
which is defined as follows (see Appendix A) 



Dij{t) = <St{t)S-{0)> 



III. DISCUSSION 



(12) 



As is well-known [18] , there can be many types of long- 
time behavior of D{t). Two of them are of particular 
interest: 



lim D{t) = L{T, H)^0 (/) 



and 



lim D{t) oc exp{—t/T) 



(11) 



(13) 
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A simple example of the system belonging to class I (so- 
called nonergodic state) has been discussed by de Gennes 
and Tinkham [19]. They considered the long-time behav- 
ior of D{t) for a superconducting thin film (of thickness 
a) with diffuse reflection on the boundaries, in a parallel 
magnetic field. When no volume defects are taken into 
account (pure limit), the system was found to exhibit a 
nonergodic behavior with field-dependent nonergodicity 
parameter L{T,H), namely 



L{T,H) 



1 - (nHayScl^o), 
{24>o/T^Ha^f, 



H « (</.o/a2) 
H > (0o/a2) 



(15) 



At the same time, the presence inside a sample of a few 
scattering centers (which are characterized by a mean- 
free path, I > a) inspires the transition of the system 
from nonergodic (class /) to ergodic (class II) state with 
the (inverse) relaxation time {vp is the Fermi velocity) 



1 / TiHa 

- = [vFa) — — 



(16) 



where 



Turning to the HTS single crystals, let us consider dy- 
namic (relaxation) and equilibrium properties of the 
magnetization versus oxygen defect concentration within 
the SG model. By analogy with the case of slightly de- 
fected thin films, considered by dc Gennes and Tinkham 
[18], we assume that up to some critical value of oxy- 
gen deficiency, 6g, HTS single crystals exhibit nonergodic 
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(phase-coherent) behavior, while for oxygen-defect con- 
centration greater than Sg, the above-mentioned coher- 
ence (within CuO plane) is destroyed and crystal un- 
dergoes the phase transition to ergodic (paracoherent) 
state where oxygen-rich superconducting "grains" are 
separated by oxygen-poor insulating boundaries so that 
there is no superconducting path through a sample. It 
is worthwhile to mention that the related problem of an- 
nealed Ising magnet on percolation clusters has been re- 
cently considered by Kaufman and Touma [20]. Using the 
renorm-group method, three phases on the correspond- 
ing phase diagram have been identified [20]: percolating 
ferromagnetic, percolating paramagnetic, and nonperco- 
lating paramagnetic. 

In view of Eqs.(10)-(13), the equilibrium magnetiza- 
tion Meq{6, T, H, Hi) is the limit 

Meg{5, T, H, Hi) = lim M{t) (17) 

t — ^OO 

= M{5,T,H,Hi)L\5,T,H). 

Here L{S,T,H) is the order parameter of the SG model 
which is defined via the correlator D{t) according to 
Eq.(13). To find the long-time (low-frequency) behav- 
ior of the correlator D{t) (and thus of the magnetization 
M(t)), we need the equation of motion for the Josephson 
pseudospins Sf{t). An approximate (vahd for 1) 
equation of motion reads (see Appendix B) [7,12-16] 

N 

s+ = pnJ2JijS+ (18) 

3 

Here Q = 2e'^R/l3h'^N is a characteristic frequency of the 
J J network with R being the resistance between grains in 
their normal state. In the so-called "mode-coupling ap- 
proximation" [21], D{t) obeys the self-consistent master 
equation (see Appendix B) 



with K{t) = {^/N)J2ij Kij{t) being a memory (feed- 
back) kernel. When there is no temporal correlations be- 
tween grains ("paracoherent state") the memory kernel 
has a "white noise" form K{t) = Kr{t) = 2^l5{t), where 
S{t) is the Dirac delta function. In this case the master 
equation results in a Dcbye-likc decay of uncorrelatcd 
paracoherent state, namely D{t) = exp^—t/r), where 
1/r = O. Such a situation is realized above some critical 
(phase-locking) temperature Tg when the coherent state 
within the J J network is destroyed completely, so that the 
order parameter L = 0. Below Tg, the situation changes 
drastically due to the superconducting correlations oc- 
curring between grains. For A^ ^ 1, the coherent part of 
the memory kernel, Kfj{t), can be approximated by the 

"current-current" correlator Kf-{t) ^ < S:^{t)Sr{0) >. 



Taking into account the equation of motion (18), the 
memory kernel below Tg can be presented in the form 
(see Appendix B) 

1 ^ 

K{t)^Kr{t) + -^Kfj{t) (20) 

ij 

= 2nS{t) + nl^i5,T,H)D{t). 

Here Q.coh{5,T,H) = (3^J{5,T,H) and the field depen- 
dence of the Josephson energy is defined as follows (see 
Appendix A) 



J{5,T,H) = Jij{5,T,H) (21) 
= JiS,T){l^^y 

In view of Eq.(13), a zero frequency (t oo) solution of 
the master Eq.(19) with the memory kernel (20) results 
in the nontrivial order parameter for the intragranular 
JJ network [7,12-16] 

The phase- locking temperature Tg{S,H), below which 
the ensemble of grains undergoes the phase transition 
into the coherent state, is defined by the equation 
L{S,Tg,H) = which, due to Eq.(22), gives rise to im- 
plicit equation, viz. Tg{5,H) = J{S,Tg, H)/kB- The 
Josephson energy depends on the temperature through 
the Ambegaokar-Baratoff relation, which near the single 
grain superconducting temperature, T^, reads J{T) « 
J(0)(1 — T/Tc). Assuming that for high magnetic fields 
(when frustration is strong enough) J{5, 0, H) ^ fc^Tc < 
J((5,0,0), we get finally Tg{S,H) « J{S,0,H)/kB. As a 
result, the order parameter L = 1 — [T/Tg{S, H)]'^ grad- 
ually changes from at T > Tg{6, iJ) to 1 at T = 0, thus 
describing the continuous phase transition. 

By analogy with the critical (phase-locking) tempera- 
ture Tg{S, H), we can introduce the critical field Hg{5, T) 
as the solution of the equation L{5, T, Hg) = 0. In view of 
the field dependence of the order parameter (see Eqs.(21) 
and (22)), the critical field reads 

Hgi5,T) = Ho^l-T/TgiS,0). (23) 

Taking into account the S dependence of the phase- 
locking temperature, Tg{S, H) ~ Tg{0, H){1 - S), Eq.(23) 
results in the following oxygen-deficiency behavior of the 
critical field 

Hg{S,T)=Ho ^5g{T,Q)-S. (24) 

Here we have introduced the critical oxygen deficiency, 
Sg{T,H), which is defined as the solution of the equa- 
tion L{6g,T,H) = and has the form 5g{T,H) = 



3 



1 — T/Tg(0, H). The physical meaning of this critical 
parameter is as follows. For S > 6g{T,H) oxygen-rich 
superconducting grains arc separated by oxygen-poor in- 
sulating boundaries so that there is no percolative path 
through the sample. Notice that within the SG model, 
Hg {6, T) in fact plays the role of the " phase boundary" 
field Hm{S,T) discussed by Osofsky et al. [3]. 

It is important to mention that the correlator D{t) fol- 
lows a simple Dcbyc-likc decay law only above Tg (S, H) , 
i.e., when the system of "grains" is in the ergodic state 
(see above). Below Tg (where the order parameter L ^ 
0), relaxation of D{t) (with Im{D{t)} = 0) can be pre- 
sented in the form [13-16] 



D{t)=L + {l-L)^t). 



(25) 



The relaxation function, $(t), is supposed to be normal- 
ized, viz. 



oo 



dt^t) = 1 



(26) 



and obeys the following boundary conditions, $(0) = 1 
and $(oc) = 0, i.e., D{0) = 1. Of course, in principle, 
one can find D{t) as a numerical solution of the mas- 
ter equation. But it seems more interesting to try and 
get some analytical results concerning the time behavior 
of D{t). It is natural then to consider a simple gener- 
alization of the Debye law in the form of the so-called 
"Kohlrausch stretched exponential" law 



^{t) = exp 



(27) 



where a{S,T,H) < 1. Substitution of Eqs.(25)-(27) into 
Eq.(19) with the kernel (20) results in an implicit equa- 
tion on the power exponent a{6, T, H) 



r( 1 + - 

a 



L 



2(1 -L)- 



(28) 



Here F is the Gamma- function. Near Tg(5,H) the ap- 
proximate solution of Eq.(28) gives 



a{6,T,H) ^1- L{S,T,H). 



(29) 



It is worthwhile to mention that for short-time limit, 
when {t/T}" <C 1, the above Kohlrausch law (27) leads 
to the non-logarithmic relaxation law for magnetization 
(see Appendix B) 

M{t) = Me,[l - 2sJog{t/T) + sllog\t/T)]. (30) 

In this approximation, the a-relaxation rate Sa{S,T,H) 
is expressed via the Kohlrausch exponent a{5, T, H) and 
the order parameter L{5, T, H) 



Sc{6,T,H)={ la. 



(31) 



In turn, Sa {S, T, H) is related to the activation energy 
U{6,T,H) as follows: Sa = ksT/U. Thus, in view of 
Eqs.(29)-(31) the S, temperature, and field dependen- 
cies of the activation energy in the JJ network are ef- 
fective via the corresponding dependencies of the order 
parameter L{d,T,H). Namely, near Tg{8,H) the acti- 
vation energy reads U{5,T,H) « kBTL{S,T,H). That 
is U{S, T, H)/2kBT « 1 - T/Tg{5, H) « 5g{T, H) - 5 Ki 
1 — H/ Hg{S,T), in at least qualitative agreement with 
what have been really observed in oxygen-depleted HTS 
single crystals [4,5] . It is interesting to notice that expres- 
sion similar to our Eq. (30) has been used by Scngupta et 
al. [22] to describe a non-logarithmic relaxation in HTS 
single crystals. 

In view of the explicit dependence of the order pa- 
rameter on the oxygen deficiency, namely L{S, T, H) w 
2{6g{T,H) — 5), Eq.(29) describes the restoration of er- 
godic state in the system under the study when oxy- 
gen deficiency 6 reaches its critical value Sg (T, H) . In- 
deed, when 5 5g{T,H), the order parameter L ^ 0, 
Kohlrausch exponent a ^ 1 (see Eq.(29)) which means 
that relaxation becomes faster (formally, according to 
Eq.(31) the logarithmic relaxation rate ^ oo) and 
follows the ordinary Debye law (see Eq.(27)). At the 
same time, the activation energy between "grains" de- 
clines (i.e., [/ — > 0) as (5 ^ Sg{T,H). To make our 
discussion more quantitative, let us consider some esti- 
mates of the model parameters. Using the experimen- 
tal results for the "phase boundary" field Hg{S,T) ob- 
tained by Osofsky et al. [3] for fixed values of S and 
T, namely Hg{5 0.06, T = 60K) w 1.5T, Hg{S = 
0.13, T = 60K) « IT, and Hg{S = 0.13, T = 70K) f« 
0.4T, Eqs.(23) and (24) allow us to get estimates for 
the phase- locking temperature Tg(6,H) and the critical 
value of the oxygen deficiency 5g{T,H). The result is: 
Tg{6 = 0,i/ = 0) « ToK, Tg{S = 0.13, iJ = 0) w 72K, 
and 5g{T = QOK, H = 0) 0.21. Finally, using the above 
results, Eq.(24) brings about the following estimate for 
the characteristic Josephson field Hq = cpo/s ^ 5T which 
gives a reasonable value of oxygen ion scattering cross 
section [3] s w 4 x 10~^^m^. On the other hand, mak- 
ing use of the above-obtained estimates we can estimate 
the value of the activation energy U{S, T, H) « 2kBT{l— 
H/Hg{S,T)). For H = IT, S = 0.06, and T = 60K, we 
get U /ks ~ 407^^, which reasonably agrees with the value 
deduced by Ossandon et al. [5] from YBa2Cu30T-s sin- 
gle crystals measurements. 

In summary, the oxygen defects induced phase transi- 
tion from the nonergodic (in nearly full-oxygenated crys- 
tals) to ergodic (in highly-deoxygcnated crystals) state in 
HTS oxygen-depleted crystals has been considered within 
the superconductive glass model. Both dynamic (relax- 
ation) and equilibrium properties of the model magne- 
tization were found to correlate quite reasonably with 
some experimental data on deoxygenated HTS. 
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APPENDIX A: 

To get Eqs.(lO) and (11) for the relaxation of mag- 
netization, wc have to calculate the relaxation function 
Using the so-called random- field approximation for 
quenched disordered systems [10,16,23], which allows to 
decouple the averaging of the "grain distribution" (repre- 
sented by the "scattering potentials" Jij) from the "car- 
riers" (or Josephson pseudospins) , i.e., assuming that 
A(ri)B(rj) = A{ri) ■ B{rj), we obtain from Eqs.(2)-(4), 
and (9) 



J2{<S+it)S-{t)S+{0)Sf{0)> (Al) 

ijkl 

X exp[i{A,j + Aki)]{xiyj - Xjyi){xkyi - xiyk) 

- exp[i(Ajj - Aki)]{xiyj - Xjyi){xkyi - xiyu) 



X < St{t)S-{t)S^{Q)S+{Q) >} + h.c. 

To proceed further, we have to calculate the 4-spin corre- 
lators appeared in the rhs of the above equation. Taking 
into account only pair correlations, namely assuming that 
(within mean-field approximation [10,16,23], when TV ^ 

1), e.g., < St{t)S-{t)S+{Q)Si{Q) >^< S+{t)Sf jO) >■ 
< S~{t)S^{0) >, we can rewrite Eq.(Al) as follows 



N 

E 

ijkl 



dH 



iAij 



dH 



(A2) 



{Du{t)D*i,{t) + Dik(.t)Djim + ^-c- 



Here we have introduced the spin-spin correlator (cf. 
Eq.(12)) 



D,jit) = <S+it)S.iO)>, 
and made use of the fact that due to Eq.(3) 

i4>o\ d 



(A3) 



{xiyj - Xjyt) exp{iAij) = ^— J ^exp(iAy ). (A4) 

The frustration field dependence of magnetization essen- 
tially depends on the choice of a random distribution 
fmiction P{fi) as well as on the type of disorder [9-15]. To 
obtain the explicit form of the field dependence of mag- 
netization given by Eq.(ll), we have assumed, for sim- 
plicity, a site-type positional disorder allowing for weak 
displacements of the grain sites from their positions of the 
original 2D lattice, i.e., within a radius d the new position 
is chosen randomly according to the normalized separa- 
ble Gaussian distribution function P(r,) = P{xi)P{yi), 
where 



P{x) = 



V2^ 



exp 



x^\ 
2d?) 



(A5) 



Using the above distribution function, we can calculate 
the configurational averages appeared in Eq.(A2). In par- 
ticular, the average value of the Josephson energy (see 
Eq.(2)) reads 



J{S, T, H) = Jij {5, T, H) = J{S, T)ex.p{iAij), (A6) 
where 

-t-oo 

7^= J dfidfjP{fi)P{fj)e^'''"/'^°'^^'''yi-''jy'^ (A7) 



1 + 772 

-"0 



2\ -1 



Here dfi = dxidyi, and Hq = (po/Trd^. Finally, taking 
into account Eqs.(A3)-(A7), we arrive to Eqs.(lO) and 
(11) for magnetization M{t). 



APPENDIX B: 



By accounting for the Kirchhoff law, lij = 0, for 
the total Josephson currents, /y = I^j + Ifj, where 
the superconducting current, 7?^, is given by Eq.(4) and 
= (fi/2eR)(d(f)ij /dt) is a normal current with R being 
the resistance between grains in their normal state, the 
approximate (valid for N ^ 1) equation of motion for 
the superconducting phase reads [10,14] 



hN d(f)i 2eJ 



N 



2eR dt 



n 



sm( 



= 



(Bl) 



Taking into account that {d/dt)e^^'^' = i{d(pi/dt)e^ 



the pseudospin representation (with 



brings 



about the approximate equation of motion (18) for 
Josephson pseudospins. 

In the so-called "mode-coupling approximation" [21], 
which is based on the Mori-like projection technique 
[23,24], the self-consistent master equation on the isother- 
mal correlation function D{t) can be constructed. Let us 
introduce the Laplace transform 



-t-oo 



(B2) 



Then the continued fraction expansion for Dij{z) leads 
to the expression [14,21] 



D,{z) = -(z 



Z + Kg{z) 



where 



1 ^ 



(B3) 



(B4) 



jk 
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Here n = 2e^kBTR/h^N is a characteristic frequency 
of the JJ network. Alternatively, using inverse Laplace 
transform, Eq.(B3) can be cast into the self-consistent 
master equation (Eq.(19)). Using the mode-coupling 
approximation scheme [14,21], the coherent part of the 
memory kernel can be represented by a set of "current- 
current" correlators 



Kt^{t) ^ < 5+(t)5r(0) > + f22< St{t)Sj{Q) > (B5) 

Since oc O (see Eq.(18)), due to a rather strong de- 
pendence of the characteristic frequency on the number 
of grains (f2 oc 1/-/V), we can restrict ourselves to a lin- 
ear approximation, Kfj{t) = < {t)Sj' {0) >, assuming 
that TV ^ 1. Taking into account the equation of motion 
(18), Kc{t) = {l/N)J2ijKfj{t) can be presented in the 
form 

Kc{t) = E E J^'^J^i ■ < (*)^r(o) > (B6) 



N 



ij kl 



= 13^9.^^(3, T,H)D{t). 

To obtain the above equation, we have used Eqs.(A3), 
(A6), and (A7) together with the decoupling approxima- 
tions discussed in Appendix A. Using Eq.(B6), we finally 
arrive at Eq.(20) for the total memory kernel below Tg. 

To get the non-logarithmic relaxation law (30) for mag- 
netization, let us rewrite Eq.(lO) taking into account 
Eqs.(25) and (27) as follows 



log(l-Z) = log 



1 — exp 



(B7) 
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